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Here in this lecture we will touch on two aspects, one the new radio meth- 
ods to observe the effects of high energy particles, and second the role that 
radio galaxies play in helping us understand high energy cosmic rays. We will 
focus here on the second topic, and just review the latest developments in the 
first. Radio measurements of the geosynchrotron radiation produced by high 
energy cosmic ray particles entering the atmosphere of the Earth as well as 
radio Cerenkov radiation coming from interactions in the Moon are another 
path; radio observations of interactions in ice at the horizon in Antarctica is 
a related attempt. Radio galaxy hot spots are prime candidates to produce 
the highest energy cosmic rays, and the corresponding shock waves in rela- 
tivistic jets emanating from nearly all black holes observed. We will review 
the arguments and the way to verify the ensuing predictions. This involves the 
definition of reliable samples of active sources, such as black holes, and galaxies 
active in star formation. The AUGER array will probably decide within the 
next few years, where the highest energy cosmic rays come from, and so frame 
the next quests, on very high energy neutrinos and perhaps other particles. 

Keywords: Cosmic rays, magnetic fields, active galactic nuclei, black holes, 
radio galaxies 



1. Introduction 

Recent years have seen a proliferation in new experimental efforts to mea- 
sure very high energy particles, both in actually constructing huge new 
arrays like AUGER in Argentina, but also new attempts to measure high 
energy particles in new ways, mostly focussing on the radio range. We can- 
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not but mention in passing the hugely successful telescopes for TeV gamma 
rays, like HESS, MAGIC, MILAGRO and Cangaroo, as well as the neutrino 
observatory IceCube for high energy neutrinos. Here we focus on yet higher 
energies, near and beyond EeV (= 10^^ eV). For the subject covered in this 
lecture, the reader is rcfcrcd to important textbooks. ^^'^ 

On the other hand, our theoretical understanding of possible production 
of high energy charged particles, neutrinos and photons is also reaching a 
measure of maturity, with many efforts concentrating on the role of shock 
waves and flows in relativistic jets, such as in radio galaxy hot spots as 
discussed in various review articles. ^"^ 

2. Radio detection methods 

It has been recognized many decades ago, that high energy particles produce 
secondary radio emission. The emission is now far along to be understood as 
geosynchrotron emission, ^"^^ or as radio Gerenkov emission. Now we have 
established efforts under way to observe, calibrate and use such emissions to 
set limits, possibly soon measure, high energy neutrinos, very high energy 
cosmic rays and also unknown particles. 

The furthest along has been the effort to use geosynchrotron emission, 
when the airshower is a directly visible radio spot in the sky. The observa- 
tion of high energy cosmic rays will be incorporated into the LOFAR array, 
in the Netherlands; the LOPES array in Karlsruhe, Germany, undertakes 
the control, and calibration of these emissions. ^^"^^ Radio Gerenkov emis- 
sion from the Moon is another effort, as is the corresponding observation 
in ice at the horizon in Antarctica. ^^"^^ There are corresponding efforts 
elsewhere, and also some tests have been done to use salt-domes. 

3. Active Galactic Nuclei 

For Galactic cosmic rays clearly supernova explosions are the prime candi- 
dates to be sources of cosmic rays; the latest results from the TeV telescopes 
strongly support this expectation and interpretation. The energy range to 
which cosmic rays derive from Galactic sources is not entirely certain, but 
the transition to cxtragalactic cosmic rays is somewhere near 3 10^^ eV, 
perhaps at slightly lower energy. 

Where ultra high energy cosmic rays come from is not certain, but the 
most promising candidates are radio galaxy hot spots and other shocks in 
relativistic jets, and gamma ray bursts, almost certainly also a phenomenon 
involving ultra-relativistic jets. 
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Here we wish to concentrate on radio galaxies, as in their case the ar- 
gument is persuasive, that they accelerate cosmic rays to extremely high 
energy. This does not necessarily imply that they are the sources for the 
events we observe. 

4. Radio galaxies 

Radio galaxy hot spots as well as many other knots in radio jets show an 
ubiquitous cutoff near 3 10^^ Hz, originally discovered in the mid 1970ies. 
Radio, infrared and optical observations strongly suggest that these hot 
spots and knots are weakly rclativistic shock waves. The basic phenomenon 
seems to be quite independent of external circumstances, and so requires a 
very simple mechanism. It has been shown that in such shocks protons can 
be accelerated,^^ subject to synchrotron and Inverse Compton losses, then 
initiate a wave-field in the turbulent plasma, and the electrons scattered 
and accelerated in the shock are then limited to emit at that synchrotron 
frequency observed, i/* ~ 3 10^^ Hz, almost independent of details. 

4.1. The maximum energy 

The maximum energy of the protons Ep^rnax from this loss limit implicated 
can then be written as 



Here we are independent of all the detailed assumptions about the inten- 
sity of the turbulence, and the exact shock speed; the dependence through 
the magnetic field on other parameters is only with the 1/7-power; typical 
magnetic field strength inferred are between 10^^ and 10^^ Gaufi. 

There is a corresponding limit from the requirement,'^^ that the Lar- 
mor motion of the particle fit into the available space; this can be written 
as Ep^^ax 

~ W^eWL%^, where Z/46 is the fiow of energy along the jet, 
some fraction of the accretion power to the black hole, in units of lO**' 
erg/s. Therefore radio-galaxies are confirmed source candidates for protons 
at energies > 10^° eV! 

4.2. Positional correlations 




(1) 



Of course it is interesting to look for associations of ultra high energy events 
and their arrival directions with known objects in the sky, whether it is the 
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supergalactic plane, some distant objects, or nearby galaxies with 
their black holes. Here we report work with loana Maris, done in the years 
up to 2004, reported on our web-page, and in many lectures. 

In order to have a statistical meaningful result we have to form com- 
plete samples of candidate sources; in some cases irregular samples will be 
unavoidable. We need active galactic nuclei, starving and active, starburst 
and normal galaxies, clusters of galaxies, with the reasoning that all such 
candidates could produce ultra high energy cosmic rays, be in shocks in jets, 
in gamma ray bursts, in accretion shocks to clusters, and hyperactive other 
stars. Recently merged black holes would be another hypothetical option, 
as then exotic particles might be produced. Most plausible would be rela- 
tivistic jets pointed at Earth, also known as flat spectrum radio sources; all 
BL Lac type sources are in this class, although not all flat spectrum radio 
sources are of BL Lac type. 

In any new search, such as to identify sources for high energy neutrinos, 
or sources for straight propagation of particles in the AUGER data, one 
should follow the same route: use complete well defined and small samples, 
and define them beforehand. With the IceCube collaboration we have just 
completed this task, and it has been published. 

We have done this task for ultra high energy cosmic rays some time ago 
with the data available publicly, using the set of accessible ultra high energy 
events with good directional accuracy above 4 IQ-'^^ eV: 80 events: AGASA 
(61), Haverah Park (6), Yakutsk (12), and Fly's Eye (1); we also used one 
event which was 38 EeV, as this was the sample used also by AGASA for 
the doublets analysis. 

This has been done before by many, from about 1960 by Ginzburg, 
including us from 1985, and in recent papers. ^^^^^ 

We find, similar to some other searches, that radio sources from the 
Condon Radio survey in positional coincidence with far infrared (IRAS) 
sources do show a highly improbable association with ultra high energy 
cosmic ray events. We can go one step further and use jet-disk symbiosis'^^ 
to predict maximum particle energy and maximum cosmic ray fiux, and 
so check whether these so identified sources can possibly produce the flux 
observed, even closely. In work with Heino Falcke, Sera Markoff, and Feng 
Yuan as well as in their own work these concepts about the physics of rela- 
tivistic jets have been tested at all levels of observability, from microquasars 
via low luminosity active galactic nuclei to more powerful sources. Current 
work is being done with Marina Kaufman. Tests include fitting the entire 
electromagnetic spectrum, and the variability. 
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Assuming that these particles can violate the GZK interaction with the 
microwave background, these sources could in fact account for the high 
flux at high particle energy, quite surprisingly.^^ Also, many of the sources 
so identified arc quite famous in radio astronomy history, sources such as 
3C120, 3C147, 4C39.25, 3C449 and the like. We did identify a speculative 
picture in which such events coming from large cosmological distances could 
be explained,^'' using particles in higher dimensions and the distortion of 
space time close to the merger of two spinning black holes. Should be here 
no such correlation once we have much more extended samples, it might be 
possible to set limits on this kind of physics. 

However, since we tried many times to get such a result we hesitate to 
assign any physical meaning to this physical picture at this time. 

5. Magnetic fields 

The magnetic fields filling the cosmos are generally too weak to influence the 
propagation of ultra high energy cosmic rays by more han a few degrees'^* 
There may be special environments where this may not be true, like the 
boundary regions around radio galaxies in a cluster of galaxies. 

However, the halos and winds of individual galaxies such as ours may 
have an appreciable influence: For starburst galaxies such as M82 the ex- 
istence of a wind has been long shown.^^''*" In our Galaxy this was finally 
recently demonstrated,^^ for NGC1808 it seems obvious in HST pictures, 
and the magnetic nature of the wind was demonstrated in the example 
of NGC4569 through radio polarization images. These last observations 
also confirm the basic ansatz**^ which we have followed in our work with 
Lauren^iu Caramete. The key point is that the bending in a magnetic field 
standard topology such as a Parker wind'''* with ~ sin^/r can be large 
as the Lorentz force is an integral in dr over 1/r, where r is the distance 
from the center, and 6 is the zenith angle in polar coordinates. Another 
more subtle point is that the turbulence in the wind is likely to be k~^, 
a saw-tooth pattern, since it is caiiscd by shock waves rimning through 
the medium produced by OB star bubbles and their subsequent supernova 
explosions. This is in contrast to the turbulence in the thick hot disk,^^ 
where observations have shown directly, that in a 3D isotropic approxima- 
tion the turbulence is of Kolmogorov nature, so is k~^^^, where k is the 
wavenumber. 

So the general concept of our Galaxy which we use, has three compo- 
nents: A cool thin disk, with lots of neutral and molecular gas, with about 
200 pc thickness, but permeated by tunnels of hot g hot disk, with 



February 5, 2008 15:54 



WSPC - Proceedings Trim Size: 9in x 6in ms 



6 

low density and high temperature with a thickness of about 4 kpc, and a 
magnetic wind, to a very first approximation akin to a Parker wind, ex- 
tending possibly quite far out, such as a few hundred kpc, or even a Mpc 
(from ram pressure arguments). 

In work with Alex Curutiu wc have calculated in a first approximation 
the scattering of ultra high energy cosmic rays in such an environment, 
assuming various possible sources, such as M87 or Cen A. We have tested 
different assumptions about the level of turbulence, and find that maximal 
turbulence alone would reproduce the homogeneous sky distribution which 
has been observed at 30 EeV. 

The predicted sky distribution (see the 2004 report on our web-page) 
consists then of long irregular stripes across the sky, as M87 is quite close to 
the Galactic North pole. One of these stripes is roughly in the same region 
as the supergalactic plane in the Southern and part of the Northern sky. 
This has been our prediction for some years now. 



5.1. Nearby sources 

In other work with Oana Ta§cau, Ralph Engcl, Heino Falcke, Ralf Ulrich, 
and Todor Stanev we have identified all available data on nearby black 
holes, and calculated their cosmic ray maximal contribution, in particle 
energy and in flux. For most sources the maximal energy of the particles is 
quite small, below or near EcV, and the flux is also low. However, for a few 
sources, such as M87, M84, Cen A, and NGC1068 the flux is of interest. At 
the highest energy only M87 competes, as already suggested many years ago 
by Ginzburg,*^ and Watson, and in a detailed physical model in Biermann 
& Strittmatter.26 

The list has also been available on our web-page since 2004. 

There is a question, however, already noted above, that a source may 

be strong, but how many of these particles make it out of the rclativistic 
bubble (visible in low frequency radio data) around the radio galaxy,''^ and 
and then on towards us? To fit the data as compiled in the PDG report*® 
the flux from M87 has to be reduced by about a factor 4, and so at 1 EeV 
Cen A is about 20 times stronger than M87, but does not itself extend 
much beyond 10^^ eV. 

Then just adding the contributions from the strongest sources, and run- 
ning them through a Monte-Carlo for propagation simulation reproduces 
the observed spectrum quite well. 
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5.1.1. Samples for testing 

Above wc have used those active black holes, for which we have data. But 
other approaches are also possible and would need to be tested. Such sample 
definitions have recently been developed in collaboration with the IceCube 
collaboration.^-'^ Here the sample selection has two key diflcrcnccs: First, 
ultra high energy cosmic rays are expected to be protons, so they interact 
with the microwave background, and so almost certainly come from nearby 
in the cosmos, about 50 Mpc or less. Second, protons are charged and so 
they may deviate from a straight line in their propagation. Here we focus 
on the sources. 

Even among the nearby sources, there are many possible candidates, 
one could consider: 

• Black holes are quite common, and so almost all galaxies have a 
massive black hole.^'^ Usually the activity of such a black hole is 
very limited, but experience demonstrates that it is almost always 
detectable in radio emission, which we interpret as the emission 
from a relativistic jet.^^ The accretion rate to power this jet could 
derive from just the wind of a neighboring red giant star. Also, 
the data suggest that black holes in the centers of galaxies always 
have a mass larger than about 10^ or 10^ solar masses;^^ we have 
argued that this minimum mass derives from the first growth of a 
dark matter star,''^'"'''* possible if in a merger of two galaxies the 
central dark matter density diverges, builds up a degenerate con- 
figuration, a dark matter star, which can then be eaten by a stellar 
black hole. Dark matter accretion has no Eddington limit, only an 
angular momentum transport limit, and that should be efficient 
during the highly disturbed situation of a merger of two galaxies. 
On the other hand, there might be some subtlety about black hole 
physics, that we are missing, and so we should obviously just check. 
The largest challenge to our physical understanding will appear if 
we find a correlation with very low mass black holes and ultra high 
energy events. At first, however, this implies just a tally of black 
holes in our cosmic neighborhood, ordered by mass. And as the 
mass of the black hole directly relates to the mass of the spheroidal 
population of stars, or the bulge, as well as the stellar central ve- 
locity dispersion, we have to start the search with these properties 
of galaxies. 

• Active black holes, such as M87, NGC315, NGC5128 (= Gen A), all 
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have rclativistic jets, and can accelerate particles. Since the power 
and also the maximum article energy scale with the radio emission, 
one needs a sample of black holes with measured compact emission, 
and so we need a sample ordered by predicted cosmic ray flux at 
energies beyond about a few EcV. This is what we attempted above. 

• If wind supernovae such as exploding Wolf Rayet stars manage to 
accelerate particles to 3 EeV, which happens to be the cutoff energy 
also for Galactic confinement, why not even further? These stars 
do not know about the Galaxy. And if some of these stars under 
unknown special circumstances produce Gamma Ray Bursts, the 
particle energy might also bo much higher. So this implies a sample 
of galaxies strong in the far infrared, where strong star formation 
galaxies, or starburst galaxies, are very prominent. If the flux in 
cosmic rays scales monotonically with the star formation rate, then 
a sample; of galaxies ordered by flux density at 60 micron would be 
the best sample to test. 

• If the propagation is delayed only slightly in its interaction with 
the intergalactic magnetic field, then also a direct connection with 
recent Gamma Ray Bursts might be worth investigating. 

• Accretion shocks to clusters of galaxies can also accelerate par- 
ticles,^^ although probably not to 100 EeV; however, we should 
actually check with observations. Most clusters, however, are very 
distant, and such a concept would qualify best for the Virgo clus- 
ter, again identifying just one most likely object in the sky, just 
like the radio galaxy M87, which is one of the dominant galaxies 
in the Virgo cluster. 

6. Predictions 

We have presented a theory to account for the entire cosmic ray spectrum 
beyond the GZK cutoff. This proposal is based on a physical and tested 
model for rclativistic jets. 

There is a Galactic magnetic wind, driven by the normal cosmic rays. 
This wind extends to some fraction of a Mpc. The existence of a Galactic 
wind in our Galaxy is now confirmed. The basic magnetic ficild topology is 
probably of Parker type. The turbulence in the wind is probably sawtooth 
pattern, i.e. and its turbulence is close to maximal. All galaxies with an 
appreciable level of star formation have such a wind, and their environment 
should look like Swiss cheese, embedded in the supergalactic sheets (work 
by L. Caramete). 
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The only contributor for cosmic ray particles beyond the GZK cutoff 
is M87, with Cen A very strongly contributing just below that character- 
istic energy, with a small contribution from NGC1068. Weaker sources are 
negligible due to their low maximum particle energy, and also due to their 
small flux. 

The arrival directions on the sky are smooth around 30 EeV, and begin 
to become patchy at higher energies, showing some characteristic stripes, 
in a very simple Parker type model for the magnetic field topology of our 
Galactic wind. 

If the arrival directions are smooth to the highest energies, then this 
source model fails, and we require Lorentz Invariance Violation, new par- 
ticles, topological defect or relic decay,^'' dark matter decay or possibly 
hints of quantum gravity. 

Should the proposal be confirmed we can develop sources such 3 C147 
as testbeds for particle physics - a CERN / Stanford / Fermilab in the sky. 

7. Acknowledgement 

P.L. Biermann would like to acknowledge Eun-Joo Ahn, Julia Becker, Venya 
Berezinsky, Geoff Bicknell, Genadi Bisnovatyi-Kogan, Sabrina Casanova, 
Mihacla Chirvasa, Alina Donea, Ralph Engel, Torstcn Ensslin, Heino Fal- 
cke, Paul Frampton, Cristina Galea, Laszlo Gcrgely, Andreas Gross, Fran- 
cis Halzen, Alejandra Kandus, Hycsung Kang, Marina Kaufman-Bernardo, 
Gopal Krishna, Phil Kronbcrg, Alex Kusenko, Norbert Langer, Hyesook 
Lee, Sera Markoff, Gustavo Medina- Tanco, Athina Mcli, Sergcj Moisccnko, 
Biman Nath, Angela Olinto, Adrian Popescu, Ray Protheroe, Giovanna 
Pugliese, Wolfgang Rhode, Sorin Roman, Gustavo Romero, Dongsu Ryu, 
Norma Sanchez, Gerd Schafer, Eun-Suk Seo, Maury Shapiro, Ramin Sina, 
Todor Stanev, Jaroslaw Stasielak, Samvel Ter-Antonyan, Valeriu Tudose, 
Ralf Ulrich, Marek Urbanik, Ana Vasile, Hector de Vega, Yiping Wang, 
Alan Watson, John Wefel, Stefan Westerhoff, Paul Wiita, Arno Witzel, 
Gaurang Yodh, Feng Yuan, Cao Zhcn, Christian Zier, now T. Kellmann, I. 
Du^an, L. Caramete, A. Curu^iu, Alina Istrate, ... 

Special support comes from the European Union Sokrates / Erasmus 
grants in collaboration with East-European Universities, with partners T. 
Zwittcr (Ljubljana, Slovenia), L. Gergely (Szeged, Hungary), M. Ostrowski 
(Cracow, Poland), K. Petrovay (Budapest, Hungary), A. Petrusel (Cluj- 
Napoca, Romania), and M.V. Rusu (Bucharest, Romania). 

Work with PLB is supported through the AUGER theory and mem- 
bership grant 05 CU 5PD1/2 via DESY/BMBF (Germany), and VIHKOS 



February 5, 2008 15:54 



WSPC - Proceedings Trim Size: 9in x 6in ms 



10 

through the FZ Karlsruhe. 



References 

1. V.S. Berezinskii, et al., Astrophysics of Cosmic Rays, North-Holland, Amster- 
dam (especially chapter IV) (1990). 

2. T. K. Gaisser Cosmic Rays and Particle Physics, Cambridge Univ. Press 

(1990) 

3. T. Stanev, High energy cosmic rays, Springer-Praxis books in astrophysics and 
astronomy. Chichester, UK: Springer, 2004. 

4. P. L. Biermann, Proc. 23rd International Conference on Cosmic Rays, in Proc. 
"Invited, Rapporteur and Highlight papers"; Eds. D.A. Leahy et al., World 
Scientific, Singapore, 1994, p. 45 

5. F. Halzen and D. Hooper, Rep. Prog. Phys. 65, 1025 (2002) 

6. T. Piran, Phys. Rep. 314, 575 (1999) 

7. G. Sigl in "Observing Ultrahigh Energy Cosmic Rays from Space and Earth" , 
Edited by Humberto Salazar, Luis Villaseor and Arnulfo Zepeda, AIP Con- 
ference Proceedings, 566, 266 - 283 (2001) 

8. B. Wiebel-Sooth and P. L. Biermann, in Landolt-Bornstein, Handbook of 
Physics, Springer Publ. Comp., p. 37 - 91, 1999 

9. H. Falcke and P. Gorham, AstropaH. Phys. 19, 477 (2003) 

10. H. Falcke, P. Gorham and R.J. Protheroe New Astron. Rev., 48, 1487 (2004) 

11. H. Falcke et al., LOPES Coll., Nature , 435, 313 (2005) 

12. W.D. Apel et al.„ AstropaH. Phys. 26, 332 (2006) astro-ph/0607495 

13. Gemmeke, H., et al., Inter.Joum. Mod. Phys. A, 21, 242 (2006) 

14. A. Haungs, LOPES Collaboration, Proceedings of ARENA 06, June 2006, 
University of Northumbria, UK, (2006a), astro-ph/0610553 

15. A. Haungs et al.. Int. Jour. Mod. Phys. A 21,182 (2006) 

16. J. R. Horandel et al., J. of Phys.: Conf. Ser., 39, 463 (2006) 

17. A. Horneffer et al. Int. J. of Mod. Phys. A, 21, 168 (2006) 

18. T. Huege and H. Falcke, Astropart. Phys. , 24, 116 (2005) 

19. T. Huege, R. Ulrich and R. Engel, Astropart. Phys. (submitted), astro- 
ph/0611742 (2006) 

20. S. Nehls et al.. Int. J. of Mod. Phys. A, 21, 187 - 191 (2006). 

21. J. Petrovic et al. Jour, of Phys.: Conf. Ser., 39, 471 (2006) 

22. S.W. Barwick et al., Phys. Rev. Lett. 96, 171101 (2006) 

23. S. W. Barwick et al., 2006, ANITA collaboration, hep-ex/0611008 

24. D. Saltzberg et al.. Intern. Journ. of Mod. Phys. A 21, 252 (2006) 

25. P. Gorham, et al., Nucl. Instr. and Meth. in Phys. Res. A 490, 476 (2002) 

26. P.L. Biermann and P. A. Strittmatter, Astrophys. J. , 322, 643 (1987) 

27. A.M. Hillas, Annual Rev. of Astron. & Astrophys. 22, 425 (1984). 

28. G. de Vaucouleurs, Vistas in Astron. 2, 1584 (1956) 

29. T. Stanev, P.L. Biermann, E. Lloyd et al. , Phys. Rev. Lett. 75, 3056 (1995) 

30. G.R. Farrar, and P.L. Biermann, Phys. Rev. Lett. 81,3579 (1998) 

31. A. Achterberg,: IceCube Collaboration and P. L. Biermann, Astropart. Phys. 
(2006, in press), astro-ph/0609534 



February 5, 2008 15:54 WSPC - Proceedings Trim Size: 9in x 6in ms 



11 

32. P. G. Tinyakov, and I.I. Tkachev, J. of Exp. and Theor. Phys. Lett, 74, 445 
(2001) 

33. N.W. Evans, F. Ferrer and S. Sarkar, Phys. Rev. D 67, 103005 (2003) 

34. Ch.B. Finley, and St. Westerhoff, Astropart. Phys. 21, 359 (2004) 

35. H. Falcke, and P.L. Biermann Astrophys. J. 293 665 (1995) 

36. Biermann, P.L., 2001 - 2006, http://www.mpifr-bonn.mpg.de/div/theory: re- 
ports 2001-2004 

37. P.L. Biermann and P.H. Frarnpton, Phys. Lett. B, 634, 125 (2006) 

38. D. Ryu, H. Kang, and P. L. Biermann, Astron. & Astroph. 335, 19 (1998) 

39. P.P. Kronberg, P.L. Biermann and F. R. Schwab, Astrophys. J. 291, 693 

(1985) 

40. T.M. Hcckman, L. Armus and G.K. Miloy, Astron. J. 93, 276 (1987) 

41. T. Wcstmcicr, C. Brns and J. Kerp, Astron. & Astroph. 432, 937 (2005) 

42. K.T. Chyzy et al. Astron. & Astroph. 447, 465 (2006) 

43. D. Breitschwerdt, J. F. McKenzie and H.J. Volk, Astron. & Astroph. 245, 

79 (1991) 

44. E. N. Parker Astrophys. J. 128, 664 (1958). 

45. S. L. Snowden et al., Astrophys. J. 485, 125 (1997) 

46. D.P. Cox and B.W. Smith, Astrophys. J. Letters 189, L105 (1974) 

47. V.L. Ginzburg and S.I. Syrovatskii, The origin of cosmic rays, Pcrgamon 
Press, Oxford (1964), Russian edition (1963). 

48. F.N, Owen, J.A. Eilek, and N.E. Kassim, Astrophys. J. 543, 611 (2000) 

49. T. Gaisser and T. Stanev, T., Phys. Rev. D 66, id. 010001 (2002) 

50. S.M. Fabcr et al., Astron. J. 114, 1771 (1997) 

51. R. Chini, E. Kreysa and P. L. Biermann, Astron. & Astroph. 219, 87 (1989) 

52. J.E. Greene, A.J. Barth and L. C. Ho, New Astron. Rev. 50, 739 (2006) 

53. F. Munyaneza and P.L. Biermann, Astron. & Astroph. 436,805 (2005) 

54. F. Munyaneza and P.L. Biermann, Astron. & Astroph. 458, L9 (2006) 

55. H. Kang, J. P. Rachen, and P.L. Biermann, Month. Not. Roy. Astr. Soc. 286, 
257 (1997) 

56. G. Amolino-Camolia and T. Piran, Phys. Rev. D 64, 036005 (2001) 

57. P. Bhattacharjee and G. Sigl, Phys. Rep., 327, 109 (2000) 



